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Abstract 


The paper is devoted to theoretical studies of the influence of cloud inhomogeneities on cloud top height (CTH) 
retrievals based on top-of-atmosphere nadir reflectance observations in the oxygen A-band. A three-demensional (3D) 
Monte Carlo code is used to simulate highly resolved spectral measurements in the oxygen A-band. These synthetic 
radiances are used as input for the retrieval code SACURA based on asymptotic radiative transfer theory and the 
independent pixel approximation. The results show that the effect of cloud inhomogeneity on the derived CTHs is small. 
While we found considerable 3D effects in the reflectance of more than 30% compared to the independent column 
approximation, the spectral dependence of the difference was small. As SACURA is mainly based on spectral ratios, the 
retrieval results are hardly affected by the large absolute deviations. In consequence, SACURA is capable to retrieve CTHs 
with an accuracy of better than 1.5km for overcast and also most partially cloudy cases. 
© 2006 Elsevier Ltd. All rights reserved. 
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1. Introduction 


The aim of this paper is to study the influence of cloud horizontal inhomogeneity (e.g., broken cloud fields) 
on the retrieval of cloud parameters using Global Ozone Monitoring Experiment (GOME) and Scanning 
Imaging Absorption Spectrometer for Atmospheric CHartographY (SCIAMACHY) (e.g., [1,3]). These 
instruments were conceived to yield global distributions of atmospheric constituents and thereby improve the 
knowledge of a variety of regional and global issues of importance for the chemistry and physics of the Earth’s 
atmosphere. In addition to trace gas columns, they allow the quantitative determination of the key 
components of the hydrological cycle, water vapor and clouds, with one single instrument. 
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To achieve the required combination of high spectral resolution and high radiometric accuracy, a large 
instantaneous field of view (IFOV) is required, resulting in a pixel size of 40 x 320 km” for GOME and 
30 x 60 km? (at minimum) for SCIAMACHY. As a result of the large IFOV, a correct treatment of cloud 
inhomogeneity is vital to exploit the full potential of GOME and SCIAMACHY. However, current remote 
sensing methods for trace gases and aerosols generally assume either cloudless sky or a homogeneous, plane- 
parallel cloud cover. In particular, cloud property retrievals explicitly or implicitly assume horizontal 
homogeneity of the cloud field within the IFOV of the instrument. On the other hand it is well known that 
clouds are highly inhomogeneous on all spatial scales, ranging from few centimeters up to the synoptic scale of 
several thousand kilometers. 

This paper is mainly devoted to the study of the influence of sub-pixel cloud inhomogeneity on cloud top 
height (CTH) retrievals based on top-of-atmosphere (TOA) nadir reflectance measurements in the oxygen 
A-band (wavelength 2 = 755—775 nm). The case of multi-layered cloudiness is excluded from our analysis. A 
review of modern techniques to retrieve CTHs is given by Rozanov and Kokhanovsky [32]. The independent 
pixel approximation (IPA) of the radiative transfer theory is commonly used for the determination of CTH. 
The IPA assumes that the reflectance of a broken cloud field can be calculated as a weighted average of 1D 
calculations for the cloudy and cloud-free parts. In particular, the IPA has been applied for the determination 
of CTH by [9,11-13,18,20], to name only a few. However, these authors just applied IPA for the problem at 
hand but they did not study its accuracy. Clearly, the applicability of the IPA for the determination of CTH 
and errors due to the usage of this approximation can be studied only with three-dimensional (3D) radiative 
transfer calculations of TOA reflectance. For our application we need spectrally resolved 3D calculations of 
the reflectance in the oxygen A-band. To our knowledge, such calculations have not been performed for 
broken cloud conditions so far. We used a well-characterized and tested 3D radiative transfer model (RTM) to 
calculate TOA reflectance in the oxygen A-band. This way we simulated satellite observations for a set of 
idealized but typical cloud scenarios and checked the applicability of the IPA for the determination of cloud 
altitudes from space. 

In the next two sections we describe the forward and inverse models used in this study. The fourth section 
discusses the uncertainty of CTH retrievals for partially cloudy scenes. 


2. The forward model 
2.1. MYSTIC 


The forward radiative transfer calculations in this paper have been performed using the 3D Monte Carlo 
code for the phYSically correct Tracing of photons In Cloudy atmospheres, MYSTIC [22,23]. MYSTIC 
is a forward Monte Carlo code which traces photons on their individual paths through a vertically and 
horizontally inhomogeneous atmosphere. The code is based on an approach similar to that described by 
Cahalan et al. [4]. Radiances are calculated using the local estimate technique (e.g., [6,21]). MYSTIC has 
been successfully validated in the intercomparison of 3D radiation codes I3RC [5]. It is operated as part of 
the libRadtran package [25] which prepares the input properties for the solver and processes the output. 
libRadtran provides several radiative transfer equation solvers and allows comparisons between those 
using identical input conditions. This way we found an agreement with the well-tested DISORT code [35] of 
better than 0.1% for various 1D cases. Excellent agreement for a 3D case very similar to this paper was also 
found by comparison with the RADUGA code, a successive orders of scattering approach for the radiative 
transfer equation, see plots in Nikolaeva et al. [29]. As for any statistical model, the Monte Carlo results are 
affected by statistical noise which decreases with the square root of the number of photons. The accuracy 
which can be reached therefore depends on the available computing power. 

For the application in this paper, we assumed a single layer cloud embedded in an atmosphere with 
vertically varying aerosol scattering and absorption characteristics (see next section). To confirm the 
performance of MYSTIC for the cases studied in this paper, we compared the results of this code with other 
radiative codes such as SCIATRAN [31] and RADUGA [29]. 
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2.2. Radiance simulations in the oxygen A-band 


To consider absorption in the oxygen A-band, a correlated-k approximation is required. For that 
purpose we use the approach of SCIATRAN which has been designed as the forward model for the 
retrieval of atmospheric constituents from SCIAMACHY [1]. Version 2.0 is used in this paper, see Rozanov 
et al. [31]. 

The SCIATRAN RTM is the improved version of the well-known GOMETRAN model which was 
originally developed to simulate solar radiation backscattered from the atmosphere and reflected from the 
Earth’s surface in the spectral range 240-800nm as measured by GOME in nadir viewing geometry. 
SCIATRAN was extended to 240-2380 nm covering the full SCIAMACHY spectral range. A new generation 
of the SCIATRAN model (version 2.0) comprises all features of SCIATRAN 1.2 and supports additionally 
radiative transfer calculations in a spherical atmosphere. The wavelength range has been extended to 175nm 
including the Schumann-Runge and Herzberg absorption bands of oxygen. The code is freely available for 
non-commercial purposes at www.iup.physik.uni-bremen.de/sciatran. The SCIATRAN 2.0 radiative transfer 
solver is based on the discrete ordinate approach similar to that described by Siewert [34]. 

For our application all relevant atmospheric processes were considered: Rayleigh scattering, aerosol 
scattering and absorption, molecular absorption, and cloud scattering. A surface albedo of 0 was assumed for 
the calculations which is a good choice over dark surfaces such as water. The model atmosphere contained 4 
aerosol layers, from 0-2, 2-10, 10-30, and 30-60km (boundary layer, tropospheric, lower and upper 
stratospheric) with different aerosol properties attributed to each layer as specified by Kneizys et al. [10]. The 
vertically integrated aerosol optical thickness t was set to 0.26 at 760 nm. A Henyey-Greenstein phase function 
pO = (1 — 9?) — 2gy + g? was used for the aerosol. Here y is the cosine of the scattering angle 0 and g is 
the asymmetry parameter assumed to be 0.6 in this work. Due to the weaker scattering of light by aerosols 
compared to clouds, assumptions on the aerosol characteristics (in particular optical thickness, single 
scattering albedo, and asymmetry parameter) do not affect the conclusions reached in this paper. 

Water clouds with different geometrical and optical thicknesses were positioned at various levels in the 
atmosphere. A constant cloud scattering phase function was assumed which was calculated using Mie theory 
at 760 nm assuming the following droplet size distribution [15]: f(a) = Aaf exp(—1.5a), where A = (3/2)’/6! is 
a normalization constant and a is the droplet radius. The refractive index of water was set to 1.33 in the Mie 
calculations. Molecular absorption was considered using a five-point correlated k-distribution [2] with a 
spectral resolution of 0.05nm, based on the HITRAN database [30]. Further details of the model used are 
given by Kokhanovsky and Rozanov [16]. 

To assure accurate results in the oxygen A-band, we compared results of MYSTIC and SCIATRAN at the 
full spectral resolution of 0.05 nm in Fig. 1 for a homogeneous cloud with optical thickness t = 20. Except for 
the statistical noise of the Monte Carlo model, good agreement is found between both codes. The systematic 
difference is very small and negligible for our application. We conclude that for 1D cases both models produce 
consistent results which is a pre-requisite for this study. 


3. The inverse model 


SACURA [15,32,33] has been developed to retrieve cloud top altitudes from TOA nadir reflectance 
observations. Results of retrievals for all SCIAMACHY data acquired globally to this date are freely available 
at www.iup.physik.uni-bremen.de. SACURA is based on the representation of the TOA reflectance R by a 
Taylor expansion with respect to the CTH h: 


R(h) = Rho) + SS ah — ho), (1) 
i=] 


where a; = R(ho)/i! and R® (ho) is the i-derivative of R at height ho. To allow inversion, the equation needs to 
be linearized [32]. We found that the function R(/) is linear over a wide range of h [16] for which reason we 
neglect the nonlinear terms in Eq. (1): 


R(h) = R(ho) + Rho yh — ho), (2) 
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Fig. 1. (Top) Reflection function in the oxygen A-band calculated by MYSTIC (open circles) and SCIATRAN (lines with points) for solar 
zenith angle 60°, nadir observation, cloud optical thickness 20, cloud geometrical thickness 0.5 km, and cloud top height 9.0 km. (Bottom) 
Relative difference between MYSTIC and SCIATRAN. 
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where R’ = dR/dh. R is measured at several wavelengths 4;(i = 1,2, ...,n) in the oxygen A-band. We therefore 
replace the scalar quantity R by the vector Rmes with components (R(1;), R(A2), ..., RUn)). Eq. (2) can then be 
written in the following vector form: 


J = ax, (3) 
where y = Raes — R(ho), @ = R (Io), and x= h— họ. Note that both measurement and model errors are 
contained in Eq. (3). A 

The solution £ of the inverse problem is obtained by minimizing the cost function ® = ||} — @x||" where |i 
is the norm in the Euclid space of the correspondent dimension. The value of £ where the function ® assumes 
its minimum can be written as 
Q, d) 


t= 
(4, a)’ 


(4) 


where (,) denotes a scalar product in Euclid space. With knowledge of the measured spectral reflection 
function Rmes and the calculated reflection function R(ho) and its derivative R’(ho) at several wavelengths, the 
value of the CTH can be found from Eq. (4) and the equality: h = £ + ho. The value of ho is set to 1.0 km, 
which is a typical value for low-level clouds. To minimize the influence of calibration errors and of a priori 
unknown atmospheric conditions, reflectances and their derivatives were normalized to their respective values 
at A* = 758nm (just outside the oxygen A-band), e.g., R(h, 2)/R(h, 2°), Rho, 2)/R ho, 2"), and (RA, A)/ 
R(h, 2*))’. The retrieved data are thus only sensitive to relative differences, rather than to absolute values. 

The reflection function R in the oxygen A-band depends not only on the CTH A but also on the cloud 
geometrical thickness / [32]. Due to multiple scattering, the pathlengths of photons are enhanced in the cloud 
and any absorber inside the cloud is more effective than outside [24]. Therefore, the value of x depends on the 
assumed cloud geometrical thickness /. Accordingly, the cost function depends on the choice of / also. To 
retrieve both cloud geometrical thickness / and cloud height h simultaneously we use the same technique as 
described above but generalized to the case of two unknown parameters [32]. 

In the next step we need a model for the spectral reflectivity R. SACURA is based on the approximate 
representation of R by the following equation [16]: 


R= Ro + Ti R.T, (6) 


where Ro gives the reflection function of the atmosphere above the cloud, R, is the cloud reflection, and 
T (i = 1,2) are the transmission from the sun to cloud top and from cloud top to the satellite, respectively. 
Approximate equations for all functions in Eq. (5) are given elsewhere [16]. The comparison of the TOA 
reflectance calculated by SCIATRAN and approximated by Eq. (5) is shown in Fig. 2 (top) for a geometrical 
thickness / = 1 km and various cloud top altitudes. As one can see, the oxygen A-band is more shallow for 
higher clouds. This can be easily understood taking into account screening of tropospheric oxygen by high 
clouds. This general feature is exploited by the CTH retrievals in this work. For the instrument spectral 
response function we assumed a Gaussian function with full width at half maximum (FWHM) 0.45 nm, which 
is close to that of SCIAMACHY. This response function is used in all calculations presented below. We found 
that the error of Eq. (5) is smaller than 5% (see Fig. 2(top)) for cloud optical thickness larger than 5. 
Therefore, Eq. (5) can be used as a basis for the retrieval procedure for the case of optically thick clouds. Thin 
clouds are out of the scope of this paper. 

The performance of SACURA for the retrieval of the CTH A and the cloud geometrical thickness / is shown 
in Fig. 3 where we present the uncertainty of the retrieved CTH A (a) and the cloud geometrical thickness / (b) 
as function of h. For this sensitivity study, the TOA reflectance has been calculated with SCIATRAN at 
l = 0.5km, /= 1km, and / = 2km assuming 100% cloud cover and various CTHs as shown in Fig. 2. All 
other conditions were the same as for the calculations shown in Fig. 2. 

We found that the CTH error is smaller than 1 km. The uncertainty of the geometrical thickness / is also 
smaller than 1 km, except for thin clouds with a geometrical thickness less than 0.5km. Such an accuracy is 
comparable to the accuracy of infrared sounders (e. g., 1.5 km as specified by [8]). The uncertainty decreases 
considerably if one of the parameters (h or /) is known in advance. In particular, one can use the correlation 
between / and cloud optical thickness t to estimate /. Such a relationship is for example given by Feigelson [7]. 
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Fig. 2. (Top) Reflection function calculated by SCIATRAN (lines) and the modified exponential approximation given by Eq. (5) 
(symbols) for a cloud with optical thickness 20 and geometrical thickness 1 km for different cloud top heights from 1 km (lowermost curve) 
to 12 km (uppermost curve) in steps of 1 km. The solar zenith angle is 60° and the observation is at nadir. (Bottom) Errors of the modified 
exponential approximation compared to the SCIATRAN simulation for cloud top heights 1, 6, and 12 km. 
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Fig. 3. Errors of the retrieved cloud top height (a) and cloud geometrical thickness (b) as functions of the cloud top height retrieved by 
SACURA for cloud geometrical thicknesses 0.5, 1.0, and 2.0 km. The solar zenith angle is 60° and the observation is at nadir. The cloud 
optical thickness is 20. 
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We obtain from their Fig. 4.2 that the extinction coefficient is typically 30 km! and get 


l= , 6 

30 (6) 
where /is in km. The value of t is easily determined from measurements outside the gaseous absorption band, 
either by a lookup-table approach or using an analytical approximation [15]. Note that Eq. (6) must be used 
with caution due to the great variability of cloud properties in space and time. Further details of SACURA are 
given by Rozanov and Kokhanovsky [32]. 


4. Retrieval of cloud parameters for broken cloud conditions 
4.1. The independent pixel approximation 


Fig. 4 shows a MYSTIC calculation of the reflectance of a simple 2D cloud. The cloud was embedded into 
the atmosphere described above. It extends horizontally from 5 to 10km in a model domain of 15km with 
periodic boundary conditions. The vertical extension was from 0.5 to 1 km, assuming a constant extinction 
coefficient 40 km™', resulting in a vertically integrated optical thickness of 20. In y-direction (perpendicular to 
the paper) the properties are constant. 

Fig. 4 illustrates the effects of 3D radiative transfer which are most pronounced for low sun: here the cloud 
sides intercept a considerably larger fraction of the radiation than a plane-parallel cloud. Part of this radiation 
is scattered into the upward direction which causes a strong peak at the sunlit side of the cloud. For high sun, 
part of the photons leaks out of the side of the clouds for which reason the reflectance is smaller than the 1D 
approximation. It is immediately clear that these effects will become more dominant for clouds which are 
“smaller” in the horizontal directions. For larger cloud elements, the reflectance approaches its 1D value 
towards the cloud center; if the cloud elements are smaller, the edge effects dominate. In fact, for cubic cloud 
elements (1 x 1 x 1 km?) we obtained differences of up to a factor of 2 between the 3D and independent pixel 
(area averaged 1D calculation, see below) results for the domain-averaged reflectivity (not shown here). 
Further discussions of these effects are given e.g. by Varnai and Marshak [37,38] and Nikolaeva et al. [29]. 

Satellite retrievals need two simplifying assumptions: (a) for technical reasons, net horizontal flux between 
individual pixels is neglected (IPA); (b) the clouds within each pixel are assumed homogeneous (plane-parallel 
approximation). While (a) is important for instruments with high spatial resolution (e.g. MODIS), (b) 
becomes more relevant for instruments with low spatial resolution like GOME and SCIAMACHY. 
Net horizontal flux may be neglected for these instruments and hence we will restrict the discussion on the 
second effect. 

Therefore, we are not so much interested in the spatial backscattered light distribution as shown in Fig. 4 
but rather in the domain-averaged reflectance. Even for the near-plane-parallel clouds used as input for Fig. 4 
we find that the domain-averaged reflectance deviates from the homogeneous value: in particular for high sun 
the reflectance is smaller than that of a homogeneous cloud due to loss of photons leaking out of the cloud 
sides and never reaching the detector. 

The spectral SCIAMACHY or GOME observations give an indication about cloud fraction within each 
pixel; however, they don’t provide any detailed information about the horizontal distribution of the clouds. A 
more detailed information may be obtained from high-resolution instruments on the same satellite platform, 
e.g. ATSR-2 on ERS-2 or AATSR on ENVISAT. Computational time constraints of operational retrievals, 
however, currently prohibit the use of more detailed information than cloud cover. Therefore we have to 
address the question how sub-pixel cloud inhomogeneity affects the TOA nadir reflectances and in particular 
the retrieval results. Assuming that we know the cloud fraction of a pixel, we test the accuracy of the IPA to 
calculate domain-average reflectances: 


R=cR.+(1—OR, (7) 


irrespectively to the horizontal distribution of the clouds. Here c is the cloud fraction which is the ratio of the 
area covered by clouds to the whole area of the ground scene studied, R, is the clear sky reflection function 
and R, is the cloud reflection function. To check the applicability of Eq. (7) we calculated R, R,, and also R, 
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Fig. 4. 3D radiative transfer calculation of the reflectance for simple 2D cloud, extending from 5-10 km in the horizontal and from 0.5 to 
1.0km in the vertical, for different solar zenith angles (sza). The position of the sun is to the left. The straight lines are the respective one- 
dimensional calculations which are basis of the independent pixel approximation. 


with MYSTIC for different cloud fractions c for a SCIAMACHY ground pixel of 60 x 30km’. As cloud 
model we adopted a random distribution of cuboids with dimensions 1 x 1 x 0.5 km? as shown in Fig. 5. The 
model resolution was 0.25km, resulting in a grid size of 240 x 120 pixels. Each cloud cuboid is randomly 
placed in the domain and cuboids were allowed to overlap. The vertically integrated optical thickness was 20 
(also for the overlapping regions). One should be aware that this is kind of a worst-case scenario: cuboid 
clouds have long been known to produce strong 3D effects, e.g. McKee and Cox [26]. Nevertheless, this 
seemingly simple cloud model is actually a good representation of a shallow cumulus cloud field. For our 
analysis we placed the clouds between 0.5 and 1 km, and between 8.5 and 9km. In particular for the latter 
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Fig. 5. Cloud model used for this study: The xy-cross section of the model domain shows a set of randomly placed cuboids, for cloud 
fractions 0.1 (top), 0.2 (middle), and 0.5 (bottom). The field contains 240 x 120 pixels of 0.25 x 0.25km?, resulting in a domain size of 
60 x 30 km? (the SCIAMACHY IFOV). Each cuboid includes 4 x 4 pixels. 


arrangement, this is a highly idealized setup. Namely, clouds at these altitudes are usually either cirrus clouds 
(which are excluded from the retrieval presented here) or high-reaching convective clouds. In order to 
demonstrate the abilities and limitations of the CTH retrieval, however, we decided to use the simplified but 
challenging cloud model described above. Fig. 6 shows the MYSTIC calculation of the reflection function for 
different cloud fractions between 0.1 and 1 for a solar zenith angle 9) = 60°. 

For the retrieval of CTH we need the reflectance of the cloudy part of the scene. Assuming that the cloud 
fraction c and the cloudless sky reflectance R, are known, we obtain from the IPA, Eq. (7): 


Re = c7! (Rmes = Ra) F Ra, (8) 


where Ries is the observed radiance. The error ¢ introduced by the IPA for Re as derived from Eq. (8) is shown 
in Fig. 7. This error is defined as e = 100(1 — R1P4/RM©), where RM is the “exact” cloudy sky reflectance, and 
RIPA is the cloudy sky reflectance approximated by Eq. (8) ( Ra is the “exact” cloudless reflectance and Rmes iS 
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Fig. 6. Reflection function in the oxygen A-band calculated by MYSTIC for different cloud fractions. 


the “exact” 3D reflectance). Fig. 7 shows this error as a function of wavelength for solar zenith angles of 0, 30, 
45, and 60°. We find that the uncertainty of the IPA for R. is strongly dependent on solar zenith angle. 
Interestingly, for this cloud geometry the uncertainty is largest for high sun which can be explained by the 
results presented in Fig. 4. Namely, for high sun and nadir observation geometry the reflectivity is reduced due 
to photons leaking out of the cloud sides. The error is largest for small cloud fractions and vanishes in the 
plane-parallel limits at c = 0 (cloudless) and c = 1 (overcast homogeneous cloud). The best agreement is found 
at solar zenith angle 60° while the differences increase again for larger angles (not shown here). The 
explanation for the small value at 60° is a cancellation of errors: the reflectance is decreased due to loss of 
photons through the cloud sides, as for high sun; this decrease is compensated by a gain due to photons 
intercepted by the illuminated cloud sides and scattered in the upward direction as shown in Fig. 4 at %9 = 75°. 
We found that errors are positive in most cases (RPA < RM), 

A very important feature to note is that in all cases the error shows only little spectral dependence. For 
retrievals that use the absolute value of the reflectance, like, e.g., for optical thickness retrievals, we expect 
large errors due to the systematic bias of the reflectance. For low sun, where the 3D reflectance is 
systematically higher than the 1D approximation (see, e.g., Fig. 4 at %9 = 75°) observations show that the 
actual reflectance may take values well outside the range of reflectivities “allowed” by the IPA [17]. SACURA, 
however, is mainly based on spectral ratios. Hence we expect the influence of the reflectance bias to be 
considerably smaller. In the next section we give a detailed analysis of the uncertainty of the SACURA- 
retrieved CTHs, by applying the retrieval to the synthetic reflectances presented in this section. 


4.2. Retrievals using the synthetic top-of-atmosphere nadir reflectance 


By applying SACURA to “exact” synthetic spectral reflectances we are able to study the uncertainty of the 
retrieved cloud products. Other than in real observations, where only the reflectances are known, here we 
know reflectances plus the real cloud properties. If the inverse problem is solved in a correct way, the 
SACURA output should coincide with the MYSTIC input. As SACURA is based on approximate radiative 
transfer theory [16] and on the IPA (see Eq. (8)), we expect some deviations. 

Fig. 8 shows the error of the SACURA~—derived CTH for a high optically thick cloud (t = 20) positioned 
between 8.5 and 9.0km as a function of cloud fraction and for solar zenith angles 0, 30, and 60°. The errors 
decrease with increasing solar zenith angle which may be explained by the smaller uncertainty of the IPA at 
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Fig. 7. Uncertainty of the independent pixel approximation for the calculation of the cloudy sky reflectance Re using Eq. (8) as compared 
to the exact result for solar zenith angles 0° (a), 30° (b), 45° (c), and 60° (d). 
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Fig. 9. Cloud top height retrieved by SACURA for a “true” cloud top height 9 km and geometrical thickness 0.5km as function of the 
cloud fraction for solar zenith angles 0°, 30°, 45°, and 60°. In contrast to Fig. 8, the cloud geometrical thickness was not prescribed using 
Eq. (6) but retrieved. 
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Fig. 10. Cloud top height retrieved by SACURA (stars) for “true” cloud top height 9km (a) and 1km (b) for a geometrical thickness 
l= 0.5km as function of the cloud fraction. Circles and squares give the values of the cloud top height, assuming a biased cloud fraction 
¢ =c+0.1. The solar zenith angle is 60° and the observation is at nadir. The “true” cloud optical thickness is 20. In the retrieval 
procedure a cloud geometrical optical thickness of 0.67 km was assumed. 


60° as shown in Fig. 7. The errors are generally smaller than 0.75 km for all solar zenith angles studied. Note 
that a constant value / = 0.67 km was assumed for the cloud geometrical thickness (see Eq. (6)) in the retrieval 
shown in Fig. 8. If the value of /is not determined independently (e.g., using Eq. (6)) but rather it is calculated 
in the retrieval procedure simultaneously with A, the errors increase as shown in Fig. 9. The different behavior 
of the accuracy with respect to the solar zenith angle shown in Fig. 9 (as compared to Fig. 8) is due to the 
compensation of the errors of the approximate equations used in SACURA [16] and the IPA (see Fig. 7). The 
conclusion from our calculations is that it is advisable to run SACURA for a fixed cloud geometrical 
thickness, e.g., as estimated from Eq. (6). Our results indicate that even if / is not precisely known, the 
uncertainty of the retrieved values of 4 is considerably reduced if the cloud geometrical thickness is available 
from independent observations. 

The cloud fraction c can be obtained using additional measurements with higher spatial resolution (e.g., the 
broad band measurements of the polarization measurements devices [36] of SCIAMACHY and GOME, or 
using ATSR-2 and AATSR placed on the same satellite platforms). In the following we show the influence of 
the uncertainties in the cloud fraction c on the CTH retrievals. Fig. 10 shows the impact of an error of 0.1 in 
the cloud fraction. The retrievals have again been performed for a fixed geometrical thickness / = 0.67 km. The 
cloud reflectance is calculated from Eq. (8). For the sensitivity study, we used the actual (“true”) cloud 
fraction as well as a biased one by +0.1. The figure illustrates that the uncertainty in the cloud fraction 
Ac = +0.1 does not lead to a considerable increase in the retrieval error for high and low clouds for the cases 
studied. Absolute errors of the order 0.1 are common in retrievals of c using, e.g., polarization measurement 
devices incorporated in GOME and SCIAMACHY [36]. 

The uncertainty of the retrieved CTH is generally larger if the value of /is determined simultaneously in the 
retrieval procedure (see Fig. 11). However, this is not related to the uncertainty in the cloud fraction but rather 
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Fig. 11. Cloud top height retrieved by SACURA (stars) for a “true” cloud top height 9 km and cloud bottom height 8.5 km as function of 
cloud fraction. Circles and squares give the values of the cloud top height (CTH) for the biased cloud fraction € = c 0.1. The solar zenith 
angle is 60° and the observation is at nadir. The “true” cloud optical thickness is 20. Results for the case A = 6km, / = 1 km are given as 
well (lower set of curves). In contrast to Fig. 10, the cloud geometrical thickness was not prescribed but retrieved. 


due to problems associated with simultaneous determination of the pair (h, /) from observations in the oxygen 
A-band alone [16]. Interestingly, inaccurate values of the cloud fraction c may lead to a compensation of 
errors of SACURA (see Fig. 11). This is related to the fact that the error in c affects mostly the absolute 
radiance but not the spectral feature used in the CTH retrievals. The cloud optical thickness can be derived 
from measurements performed outside the O) A-band. We expect much larger errors in the retrieved values of 
t compared to h for broken cloud fields. Fig. 12 shows that this is actually the case. These differences are best 
explained by the fact that the retrieval of t depends on absolute values of the reflection function while for the 
CTH determination we use ratios of reflectances and not absolute values of R. 


5. Conclusions 


Modern satellite trace gas retrieval algorithms must be improved with respect to the incorporation of 
detailed information on clouds present during the measurements. The restriction to completely cloudless pixels 
reduces the information content of measurements considerably. The main cloud parameters needed, e.g., for 
the determination of the total ozone, are the cloud fraction, the CTH, and the cloud optical thickness. These 
parameters can be determined using cloud retrieval algorithms. An important requirement for such algorithms 
for the application to studies of atmospheric chemical composition are their accuracy and speed. SACURA is 
based on asymptotic analytical equations, and, therefore, it satisfies the speed constrains. Its accuracy is also 
satisfactory for problems related to the trace gas vertical column determination. We summarize that 
SACURA is capable to retrieve cloud top altitudes not only for overcast scenes (see Fig. 3) as it was stated 
before [32] but also for partially cloudy scenes, if the cloud fraction c is obtained from independent sources, 
e.g., from the polarization measurement devices of GOME and SCIAMACHY [36]. 
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Fig. 12. Cloud optical thickness retrieved by SACURA (stars) for a “true” cloud optical thickness 20. Circles and squares give the values 
of the cloud top height assuming a modified cloud fraction c = c + 0.1 in the retrieval. The solar zenith angle is 60° and the observation is 
performed at nadir geometry. 


In addition to the standard version of SACURA [32], the cloud reflectance is determined using the IPA 
based on the cloud fraction and climatological values of the clear sky reflectance R,.. Ra may also be obtained 
from minimum values of R, for a given pixel or area, see e.g., by Koelemeijer et al. [14]. This is of special 
importance for land surfaces, where the unknown surface albedo (assumed to be zero in this study) may 
further bias results. Accurate values of R, are not needed for values of c close to one. In this case the 
reflectance is dominated by the cloudy portion of the scene. 

Our studies showed that 3D radiative transfer may have a considerable effect on the absolute values of the 
reflectance. However, the deviations between exact 3D theory and the IPA were relatively constant over the 
spectral range of the oxygen A-band. Hence, the retrieved optical thickness (based on absolute reflectance 
values) is strongly affected by 3D radiative effects while the impact on the determination of CTH (based only 
on spectral ratios) is only slightly sensitive to 3D effects. 

The error of the cloud top altitude determination using SACURA for the overcast clouds was smaller than 
l km, if both cloud geometrical thickness and CTH are determined in the retrieval procedure. The error 
depends on the cloud geometrical thickness and is larger for thinner clouds (see Fig. 3). We also found that the 
error is highly sensitive to a priori assumptions on the cloud geometrical thickness and also on the solar zenith 
angle and the cloud fraction (compare Figs. 8 and 9). It is considerably reduced, if the value of the cloud 
geometrical thickness is obtained from independent measurements. This indicates a possibility to retrieve the 
cloud geometrical thickness and also cloud bottom height from reflectance measurements in the oxygen 
A-band, e.g., if the CTH is determined by a lidar [19]. In most broken cloud cases the error of SACURA was 
smaller than 1.5km which is close to typical errors of IR measurements as applied to the determination of 
cloud altitudes from space. For instance, Frey et al. [8] reported that the CO, slicing method retrieved cloud 
heights to within 0.5 and 1.5 km of the lidar heights in 32% and 64% of the cases, respectively. Naud et al. [28] 
demonstrated that MODIS derived CTH altitudes differ from that of a ground-based radar by 1-2 km. Minnis 
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et al. [27] have shown that IR measurements are capable to determine altitudes of clouds with the error less 
than 0.5 km for thick extended clouds and to within 1.0km for thin clouds. 

We believe that the results obtained here and also the Monte-Carlo technique used are of interest for several 
scientific groups conceiving cloud altitude determination algorithms using backscattered radiation in gaseous 
absorption bands and generally for those interested in the development of accurate procedures to determine 
trace gas concentrations from orbiting and geostationary satellites. 
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